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1. Introduction and summary 

The Dyson-Schwinger equation (DSE) is a powerful tool to investigate with 
the analytic procedure the nonperturbative structure of field theories, such 
as the phase structure of gauge theories. However, the full DSEs are cou- 
pled integral cqutions for several unknown functions, thus are hard to be 
solved without introducing appropriate approximations. We usually adopt 
the step-by-step approach to this problem, firstly approximate the integra- 
tion kernel by the tree, or, ladder kernel, next use the improved ladder 
one, etc. Advantage of the DSE analysis lies in the possibility of such an 
systematic improvement through the analytic investigation. 

Analyses of the DSE have proven to be successful in studying the phase 
structure of vacuum gauge theories [1-3]. In the Landau gauge DSE with 
the ladder kernel for the fermion mass function in the vacuum QED, the 
fermion wave function renormalization constant is guaranteed to be unity 
[1], satisfying the Ward identity. Thus irrespective of the problem of the 
ladder approximation, the results obtained would be gauge invariant. 

Same analyses have been carried out in the finite temperature/density 
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case with the ladder kernel [4-6], and with the hard-thermal-loop (HTL) 
resummed improved ladder kernel [7] . Results of Ref . [7] show that at finite 
temperature/density it is important to correctly analyse the physical mass 
function the mass function of the "unstable" quasiparticle in thermal 
field theories, and also to correctly take the dominant thermal effect into 
the interaction kernel. 

All the preceding analyses [4-7], however, suffer from the serious problem 
coming from the ladder approximation of the interaction kernel. Although 
in the vacuum case, despite the use of ladder kernel, in the analysis in 
the Landau gauge the Ward identity is guaranteed to be satisfied, at finite 
temperature /density there is no such guarantee. In fact, even in the Landau 
gauge the fermion wave function renormalization constant largely deviates 
from unity [7,8], being not even real. At finite temperature/density the 
results obtained from the ladder DSE explicitly violate the Ward identity, 
thus depend on the gauge, their physical meaning being obscure. 

In this paper, we worked out, in the analysis of the HTL resummed 
improved ladder DS equation for the fermion mass function in thermal 
QED, the procedure to get the gauge invariant solution in the sense it 
satisfies the Ward identity, and investigated the properties of the "gauge 
invariant" solution. 

Results of the present analysis are summarized as follows: 

(1) We can determine the solution that satisfies the Ward identity, 
namely the fermion wave function renormalization constant being almost 
equal to unity. To get such a solution it is essential that the gauge parameter 
£ depends on the momentum of the gauge boson. 

(2) The chiral phase transition in the massless thermal QED is confirmed 
to occur through the second order transition. 

(3) Two critical exponents v and 77 are consistent with constant within 
the range of values of temperatures and couplingconstants under analysis: 
v = 0.395,7? = 0.518. 

(4) The effect of thermal fluctuation on the chiral symmetry breaking 
and/or restoration is smaller than that expected in the previous analysis in 
the Landau gauge [7]. 

2. The Dyson-Schwinger equation for the fermion 
self-energy function Ejj 

The fermion self-energy function S« appearing in the fermion propagator 
S R (P) = lf+ie 1 °-i: R (P)}- 1 (1) 
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can be decomposed at finite temperature and/or density as 

E fl (P) = (1 - A(P)) Pl Y - B(P)j° + C(P) (2) 

with A(P), B{P) and C(P) being the three scalar invariants to be deter- 
mined. In the present analysis, we use the HTL resummed form *G A " y for 
the gauge boson propagator G^ v 



t fip,v I is\ _ \ aplv j \ fi^ v — 

1 } ~ *n T - K 2 - iek Q *n L - K 2 - iek K 2 + iek Q 



(3) 



where *Hl/t 1S the HTL resummed longitudinal/transverse photon self- 
energy function [9], and A^ v , B^ v and D^ u are the projection tensors [10], 

A v» = g p» _ B ^ _ ^ BP ,v = _k^K v /K 2 , D" v = Kf*K v /K 2 , (4) 



where K = [k, fcok), k = v k 2 and k = k/fc denotes the unit three vector 
along k. The parameter £ appearing in the term proportional to the projec- 
tion tensor D pv represents the gauge-fixing parameter (£ = in the Landau 
gauge). This gauge term plays an important role in the present analysis. 

The vertex function is approximated by the tree (point) vertex. With the 
instantaneous exchange approximation for the longitudinal photon mode, 
we get the DSEs for the three invariant functions A(P), B{P) and C(P) 



-P 2 [1-A(P)] 



[ d 4 K 



{1 + 2n B (p - fco)}Im[ *G P °(P - K)\ x 



{K a P p + K p P a - p (K a g p0 + K p g a0 ) - k {P a g p0 + P p g a o) + pkzg ap 

A{K) 



+2p Q k g r7 o9po} 



[ho + B(K) + ^e] 2 - A(K) 2 k 2 - C(K) 2 
fco + B{K) 



+ {Pa9p0 



+P p g a0 2p a g a0 g p0 } ^ + b{k) + _ A(K)2k2 _ C{K)2 _ 

+{1 - 2n F (k )} *G P °(P - K)lm[{K„P p + K p P a - po{K a9p0 + K p g M ) 

-k (P a g pQ + P p g a o) +pkzg ap + 2p Q k a g a0 g pa } x 
A(K) 



[ko + B(K) + ^e] 2 - A{Kfk 2 - C{K) 2 
ox k " + B W 



+ {Pa9p0 + Pp9a0 



[k + B(K) + ^e] 2 - A{K) 2 k 2 - C(Kf 



(5) 



February 2, 2008 8:57 



WSPC - Proceedings Trim Size: 9in x 6in proc 



4 



-B(P) 



d*K 



(2tt)< 

{K a g P o + K p g a0 - 2k gaog P o} 
~{2g P o2g a o - gap} 



{1 + 2 B ( Po - fc )}Im[ *G%(P - K)] x 
A{K) 



[ko + B(K) + ^e] 2 - A(K) 2 k 2 - C(K) 2 
k + B{K) 



+{1 - 2n F (fc )} *G%(P - K)Im 

x{K a g P o + K p g a0 - 2k g a0 g P o} 
k + B(K) 



[fc + B(K) + ie} 2 - A(K) 2 k 2 - C{K) 2 

A(K) 



[k Q + B(K) + ie] 2 - A{K) 2 k 2 - C{Kf 



+ 



[ko + B(K) + ie] 2 - A{K) 2 k 2 - C(K)' 



-{2g P o2g a o- g<r P } 



(6) 



C(P) 



7 



d*K 
(2tt)4 
C(K) 



g ap {l + 2 B ( P o - fc )}Im[ *G%(P - K)] x 
+ {l-2n F (fc )} x 



[ko + B(K) + ie] 2 - A{K) 2 k 2 - C{K) 2 
"G^(P-K)Im^ C[h] 



.[ko + B{K) + ie] 2 -A{K) 2 k 2 -C{K) 2 \\ ' (?) 

The function A(P) is nothing but the inverse of the fermion wave func- 
tion renormalization constant Z2, thus must be unity in order to satisfy 
the Ward identity in the ladder DSE analysis, where the vertex function 
receives no renormalization effect, Z\ — l. 

We must solve the above DSEs and get the solution satisfying the Ward 
identity Z 2 = Z\{= 1), where Z 2 = A(P)^ 1 . The procedure to get the 
"gauge invariant" solution is as follows; 

(1) Assume the nonlinear gauge such that the gauge parameter £ to be 
a function of the photon momentum K = (ko, k), and parametrize £ as 

£(fco,k) = ^2U n H m (ko)L n (k), (8) 

where £ mrl are unknown parameters to be determined, H m the Hermite 
functions and L n the Laguerre functions. 

(2) In solving DSEs iteratively, impose the condition A(P) = 1 by con- 
straint for the input-functions at each step of the iteration. 

(3) Determine £„ m so as to minimize \A(P)— 1\ 2 for the out-put functions 
and find the solutions for B(P) and C(P). 

3. Results of the analysis "gauge invariant" solution 

Here we present the results obtained by allowing the gauge parameter £ to 
be a complex value. Number of parameters £ m „ to minimize |A(P) — 1| 2 
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is 2 x 5 x 2 = 20 (i.e., m = 1 ~ 4 and n — 1,2). All the quantities with 
the mass dimension are evaluated in the unit of A, the cut-off parameter 
introduced as usual to regularize the DSEs. 

Now we present the solution consistent with the Ward identity, i.e., the 
"auge invariant" solution. Firstly in Fig.l we show KeA(P). For comparison, 
results in the constant £ analyses are also shown in the same figure. 



a=4.0 



^=-0.05 

^=o.o • 

^=0.05 - 



"'0.12 0.125 0.13 0.135 0.14 0.145 0.15 0.155 

T/A 

Fig. 1. Comparison of the renormalization constant Re[ J 4] at the coupling constant 
a = 4.0 evaluated at po = 0, p = 0.1, see, text. 

Next let us study the property of the phase transition. Fig. 2(a) shows 
the real part of the fermion mass ReM(P), M{P) = C(P)/A(P), obtained 
from the "gauge invariant" solution, as a function of the temperature T . 
The mass is evaluated at po = 0, p = 0.1, to be consistent with the standard 
prescription to define the mass in the static limit, p = 0, p — > 0. 



1.4 
1.2 . 



0.8 
0.6 

n A 




Fig. 2. (a)Temperature-dependcncc of the fermion mass Re[M(P)] for various values 
of the coupling a. (b) Coupling constant-dependence of Rc[M(P)] for various values of 
the temperature T. Both are evaluated at po = and p = 0.1. As for the various curves, 
see text. 



The curves in the figure show the best-fit curves in determining the 
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critical temperature T c and the critical exponent v, by fitting, at each cou- 
pling constant a and near the critical temperature T c , the temperature- 
dependent data of ReM(P) to the functional form 

ReM(P) = C T (T C - Tf. (9) 

Also shown in Fig. 2(b) is the ReM(P) obtained from the "gauge in- 
variant" solution, as a function of the coupling constant a. The mass is 
evaluated at po = 0, p = 0.1 as above. The curves in the figure show the 
best-fit curves in determining the critical coupling a c and the critical ex- 
ponent -q, by fitting, at each temperature T and near the critical coupling 
a c , the coupling-dependent data of ReM(P) to the functional form 

ReM(P) = C a (a — a c ) n . (10) 

The determined critical exponents are given in Table 1. 

Table 1. Critical exponent v for various 
values of the coupling constant a and critical 
exponent r) for various values of the temper- 



ature T 






a 


V 


T 


V 


3.5 


0.42800 


0.115 


0.54718 


4.0 


0.38126 


0.120 


0.57872 


4.5 


0.36420 


0.125 


0.51430 


5.0 


0.40579 


0.130 


0.46153 



The averaged value of v over the various coupling is < fi >= 0.395, 
which fits to all the data ReM(P) in Fig. 2(a) irrespective of the value of 
the coupling constant. The averaged value of r\ over various temperatures 
is < r\ >= 0.518, which fits to all the data ReM(P) in Fig. 2(b) irrespective 
of the value of the temperature. 

The phase boundary curve in the (T, a)-plane thus determined shows 
that the region of the symmetry broken phase shrinks to the low- 
temperature-strong- coupling side compared with that of the Landau gauge. 
This fact means that the effect of thermal fluctuation on the chiral sym- 
metry breaking/restoration is smaller than that expected in the previous 
analysis in the Landau gauge [7]. 

4. Discussion and comments 

Results presented in the present paper are preliminary, because of the rough 
analysis of the data processing. We are now refining the data analysis and 
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soon get the results of the thorough reanalysis. Though the main conclusion 
will not be altered, several remarks should be added. 

(1) Present analysis was performed by allowing the gauge parameter £ 
to be a complex value. Such a choice of gauge may correspond to studying 
the non-hermitc dynamics, thus may cause some troubles. What happens 
if we restrict the gauge parameter to the real value? We are studying this 
case, finding a remarkable result: In both cases results completely agree, 
thus getting a solution totally independent of the choice of gauges. 

(2) In the present analysis, the consistency of the solution with the Ward 
identity is respected only by imposing the condition A(P) m 1. Needless to 
say, in solving the (improved) ladder Dyson-Schwinger equation, there are 
no solutions totally consistent with the Ward identity. Despite that fact, 
following point should be closely examined: At least around or in the static 
limit, po = 0, p — > 0, where we calculated (defined) the mass, each invariant 
function A(P), B{P) or C(P) should not have big momentum dependence. 
This condition may be important in connection with the consistency of the 
obtained solution with the gauge invariance. Result of the present analysis 
shows that at least B(P) and C(P) satisfy this condition. 
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